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Abstract 


Nonlinear effects of parametric and of heating type , produced in a plasma 
under the action of an electric field E = E 0 e ,u,t , are considered in this work in 
connection with the so called Tethered Magnetospheric Cloud (TMC) accom- 
panying the Tether Satellite System (TSS). The theoretical results presented 
below, particularly by some numerical calculations, show that these phenom- 
ena should appear in the ionosphere at high altitudes Z > (150 — 200) km, 
particularly, at Z ~ 300 km of the TSS system orbit. Therefore, it is 
of a special interest to search these phenomena by such a unique experi- 
ment as the forthcoming first TSS-I and by the future, perhaps modified 
TSS missions. Because of the parametric decay instability , new branches 
of waves may be excited both around the electron and ion Lengmuir fre- 
quencies ujq — 27t/ 0 and fl 0 = 2ttF 0 under the influence of high frequency 
(HF), (/ < 10 6 to few 10 6 Hz ), strong (|E| ~ V/m to tens Vfm) electric 
waves. The heating of all the kinds of particles is growing up very quickly in 
the ionosphere with altitude in the extra low and very low (ELF and VLF) 
frequency ranges F a (1 to 10 4 ) Hz discussed below. The temperatures 
(energies), for example, of the electrons accelerated by the electric field be- 
come larger than the ionization potential in this frequency range already at 
altitudes Z > (150 — 200) km when the amplitude of the electric field is 
l-Eol (1—2) mV/m. The sources of these electric field may be in the TSS-I 
mission the so called Phantom Loop (PL) - the Tethered Electrodynamic Tail 
(TET), and different kind of e.m. oscillations produced by different kind of 
instabilities in the TMC plasma. The growth rates of these instabilities will 
become very high in the TSS surrounding magnetoplasma. However, in the 
future TSS missions special artificial sources (generators) of electric fields 
should be used for these investigations. 
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Section I. Introduction 


The nonlinear behavior of a magnetoplasma, let us say, of the ionosphere, 
under the action of an outer source of an electric field E = E 0 e xut was studied 
by many authors. This problem has a long history. It began by the papers 
[1-4] (see Bonch Bruevich (1932), Tellegen (1933), Bailey and Martyn (1934), 
Forsterling (1935)) and by many others in (1933-38) - after the discovery of 
the crossmodulation of radiowaves of different frequencies, reflected from the 
ionosphere, due to the nonlinear interaction of these waves. Shortly after, it 
was clear that the mechanism of this effect is created by the increasing of the 
velocities of the electrons v e and of the other constituents of the plasma under 
the influence of the electric field. Namely, the velocities of all the constituent 
particles of the ionosphere and the collision frequencies v between them, and 
their temperatures T e ^ >n i.e. the conductivity of the plasma become functions 
of the amplitude Eo and of the angular frequency u> of the electric field E. I.e. 
v = v(EoiU>), v = v e , lyn (E 0 ,uj,v(E 0 ,Lo)) and T = T eyl , n {E 0 ,v(Eo,u)), where 
the indexes e, i, n denote, respectively, electrons, ions, and neutral parti- 
cles. Thus, all the equations which determine these values and describe their 
behavior in the magnetoplasma become nonlinear. To learn theoretically dif- 
ferent aspects of this problem, the selfconsistent solution of the adequate 
systems of equations must be solved. 

It is selfevident that these heating type phenomena , occurring in a colli- 
sional plasma are acting in large plasma regions. The influence of the electric 
field covers regions with linear scales much larger than the mean free paths 
of the particles. These effects were described in dozens of papers and in some 
monographs. Many references of these works can be found in the compre- 
hensive monograph in this field by Gurevich (1978, [5]). Earlier references 
are cited in the Russian book by Alpert (1947, [6]). 

Another kind of nonlinear phenomena are the so called parametric wave 
decay effects. They were discovered about three and more decades after the 
discovery of the heating effects (see Silin 1965 [7], Du Bois and Goldman 1967 
[8]). These phenomena are local in space and should not play a large role 
in the heating effects. The general theory of these phenomena is described 
in the comprehensive monograph by Silin 1973 [9]. These phenomena were 
studied just a little in connection with the ionosphere. The known theoretical 
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results are used below, by considering some of these effects in the Tethered 
Magnetospheric Cloud [10]. 

However, the theory of the heating effects, known from the literature 
(see, for example, [5]), is working in a limited frequency band, namely when 
u> 2 Jf + v e v in (u>L = 2 nFi, F l is the lower hybrid frequency), i.e. in the 
ionosphere at frequencies F > (4 - 5)10 4 Hz. Additionally in these studies, 
the influence of the velocities of the ions and of all neutral particles, and also 
the influence of the collision frequencies v tn between the ions and neutral 
particles, etc. were not taken into account in detail. 

In this work, results of theoretical study of the heating of a magneto- 
plasma in all the frequency range are given, taking into account all kinds 
of collision frequencies and also their temperature dependencies, i.e. their 
dependence on the electric field E = E 0 e ,we . Some results of numerical cal- 
culations of the temperatures T e , T, and T n of the electrons, ions and neutral 
particles are presented here for the ionosphere in altitude, frequency and an- 
gle 0 dependencies in the regions Z = (100 — 1000)fcm, F — (1 — 10 4 )Hz 
(0 is the angle between the wave’s electric field E 0 and the geomagnetic 
field Ho). Some of these calculations are based on the self consistent solution 
of two systems of equations in the hydrodynamic (microscopic) approxima- 
tion. One of these systems (see below (30)) consists of three vector (nine 
scalar) Lorentz equations of the derivatives of the velocities V e , V,, V n ol 
the particles. It is supposed that the magnetoplasma consists of one kind of 

neutral particles. The second system of three equations (see (31)) is of the 

dT dT 1 dT 

derivatives of the temperatures 

The general formulas of the self consistent solution are very complicated 
and bulky: they are absolutely immense. Only numerical results should be 
used to learn the T e ^ n and V e ,;,n dependencies in this case. However, in the 
v = const approximation , i.e. when v e { = v e i , o , v e n = t'en,o 5 ^tn = 
and the electric field dependence of v is not taken into account, the ana- 
lytical formulas obtained are sufficiently visible. These formulas were also 
used for some calculations. 1 It was important to evaluate the applicability of 
the formulas recommended in the literature, obtained for v — const. For this 

H am very grateful to Dr. Bob Estes who did all the calculations by computer. 
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purpose and also for the full study of this problem, the velocities and temper- 
atures are calculated, in some cases, both, with the full systems of equations 
(3 and 3) of (30) and (31), and with the shortened systems of equations (3 
and 2). In the last case, the neutral particles temperature T n = T n o = const. 

The dependencies of the temperatures on time, namely, the process of 
their transition to the stationary state is also illustrated by some examples. 
The role of the neutral particles is important in this process. The system of 
equations of the temperatures . has not a stationary solution. Math- 

ematically it is clear: the determinant of this system A = 0. Therefore T n is 
growing up in time. It means that the temperature of the neutral particles 
becomes even a source of the heating of the electrons and ions. From the 
very beginning of the action of the electric field, the source of the heating 
of the neutral particles is the energy transferred to them by the collisions 
with the electrons and ions. The values ( T e — T n ) and ( Ti — T n ) are positive. 
However, growing up in time, T n becomes equal to T t . From that moment, 
the degree of the growth of the temperatures becomes larger and larger, and 
the temperature of the electrons becomes larger than it would be without the 
influence of the neutral particles. The process as a whole has not a stationary 
solution. To stop the growth of the temperatures, in addition to collisions, 
other sources of loss of the energy accompanying the heating of the magne- 
toplasma must be taken into account, (see below Section III. 2). 


Section II. Nonlinear parametric phenomena 


The mechanism, called parametric decay instability in the strong pump 
waves , E p cos co p t can excite in a magnetoplasma resonance branches of waves 
and oscillations in the extra low (ELF, 0 < < ft//), very low (VLF, ft// < 

u> < low (LF, ujl < u) < u>//), and high (HF, lj > u >//) frequencies. These 
effects become active when the energy density of the electric field is much 
larger than the pressure of the electrons, namely: 


N-kT- 


«r- 


It means that in the ionosphere, at the altitude Z 
electron density and temperature are equal to N e ; 


( 1 ) 

300 km, where the 
1.8 • 10 6 cm 3 , T e ~ 
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1.5 • 10 3 A'°, the electric field should be rather strong , namely E p > (60 - 
70) V/m. However, it is shown below that in our case, close to the resonance 
region lo p ~ u>q, namely when 

0 < w p - y/ul + < \k 2 D 2 sf4 + ft 2 0 , (2) 

the amplitude Eq of the electric field becomes much smaller: E p is a few 
Volts / m (see Section II.3). In (1 - 2) k = ^ is the wave number, A is the 

wavelength of the parametric excited oscillations, D = ( 4 jrjVe 2 ) * s 

bye length, e and m are the charge and mass of the electrons, Ho is the ion 

Langmuir frequency, and the Boltzman’s constant k — 1.38 -10 erg -deg 


Let us note here that when we are looking for the parametric effects in 
the ionosphere or in the magnetosphere under the action of artificial sources 
E p cosu p t, located on rockets or satellites, or excited in the surrounding 
region of the magnetoplasma, the power of these sources should not be very 
large. Let us note that the growth rates 7 eo and 7 e H of excitation by beams of 
electrons of natural oscillations in the HF frequency band by the Cherenkov 
and cyclotron resonances 


7e0 

7 eH 


N b (V b \ 2 2 „ 

TT — w O COS fc> , 

N. \v e J 


4 N e v e 


( 3 ) 


are very large. By discussing this problem for the XSS mission, the values 
of the growth rates become even comparable with u>o (Alpert 1989, 1990 [10]). 


In the theory of the parametric nonlinear effects, a characteristic param- 
eter appears 


eEp 2 tt 

mio'* A 


(see Silin 1973 [9]). In the ionosphere at Z ~ 300 km 


Pv ~ “7 ("y) E p , t o p - .su>o, s > 0 
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To give a general presentation of the parametric effects by considering 
them at the altitude Z = 300 km , some formulas for an isotropic plasma can 
be used for simplicity, because Uq u?//. 


II. 1 Dispersion equations, strong electric field 
E — E p cosu; p t, uj p ujq , iv’H , Z = 300 km. 

The dispersion equation of longitudinal oscillations of a magnetoplasma, 
excited under the action of a strong electric field of the angular frequency u> p , 
which is much larger than the electronic angular Langmuir lj 0 and gyro u>// 
frequencies, not taking into account the thermal velocities of the particles oi 
the plasma, is the following: 

rw x 8 [, (Wq + Oq) COS 2 0" 

U(U!) - 


U 


1 - 


ur 


1 - 



1 - 


U 1 


+ 


+ 


u; 6 sin 2 0(^ + ^)(l-^^) + 

i — </o(/h>)] fi- 


le) 



= 0 


In (6), in addition to the notations given earlier, J 0 is the Bessel function, 
tin is the angular ion gyrofrequency, and p v - see (4), (5). 


The solution of (6), estimates, in the absence of the electric field, when 
E 0 = 0 and J${p) = 1, the wellknown three resonance branches of the plasma 
waves 


0 = 0: 


0 = 


7T 

2 


w, 


Taking into account that 
frequencies are 


2 2 2 2 ^ 

: i Itf, w 3 = H) + “0 

(l + 

: 0 , u 2 2 =u 2 l = (u; h SI h )± 

: w u = { u o + ^0 + + &h) • 

^ | C 1 the lower and upper hybrid 




( 8 ) 
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In the presence of the electric field, the solution of (8) gives four resonance 
branches 


0 


2 r\2 2 q 2 l 1 J ° 2 M ~ O 2 [l — Jn(p ) 

= 0 : — H//, u; 2fP — — 2 t n2 0 r o\Pp) 


loq + ng 

1^3 — UjJj , — Cl^q + — ^0 


( 9 ) 


and 


7 r 


0 = - : w* = 0, cj 2 = 0, 


UJ 


2 _ 2 , ^0^ 

3 = U.„-U, t + ^ 


( 10 ) 


— ufy = UOq + fto + U H + ^// + U \p — + W H + W Lp ■ 


In (6) - (10) ft//, u>l, w// are the ion gyro, low hybrid and electron gyro fre- 
quencies, ft 0 and loo are the ion and electron Lengmuir frequencies, u> v is the 
hybrid frequency and in the formula of is used the inequality ft// fto- 


Approximate theoretical dependencies of these branches of wave on the 
angle 0 P between the wave vector k and the geomagnetic field Ho at the 
altitude Z = 300 km are given for illustration on Fig. 1 both for E 0 = 0 and 
E 0 ^ 0. The characteristic values of the frequencies used in these calculations 

are 

o; 0 = 7.6 • 10 7 , ft o = 7.9-10 5 , = 7.7 ■ 10 6 , ft// = 4.14 • 10 2 

(see the Table in Section IV), and p p = 1, 0.1 and 5 • 10 4 . For excitation of 
the new parametric branch of oscillations u> 2 lP — (ftp — y 0), the amplitude Eq 
of the electric field should be about a few hundreds Volts jm when w 2 > cUq. 


By increasing p p , the frequency of ft p is also increasing. For example, 
when p p = 1, ftp ~ 0.51 ft 0 and approaches quickly the ion Langmuir fre- 
quency fto- By decreasing p p , namely when p p <C 1, Jq{Pp) — (•*■ ~ *r) an< ^ 

n, = n. [i - 4(o,)] ,/J = -^ = ? (f ) (u) 
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(see (5)). The frequency Q p becomes small and can reach ELF frequencies 
u <,<C 0// in the bound of the Alfven waves (see (5)). The amplitude of 
the electric field needed for the excitation of these very low frequency waves 
remains large when u p » Uq. However, by approaching the resonance region 
u> p ~ w 0 (see below II. 3) the electric field becomes much smaller, namely, 
E o a few Volts/m. 

The new branch of oscillations aJ 2lP = (fi p — ► 0) is of special interest. 
Taking into account the thermal velocities of the charged particles v e and t’,-, 
i.e. in the kinetic approximation , the equation of this resonance branch is 
the following: 

U>2 p = COS 2 0 




K (...) + 3k 


’-V 


(12) 


where 


K(...) = 
COS 0 


nsii-Jo a (<>,)](! + **p;r 


-1 ’ 


fi# + fig sin 2 0 [1 - J 2 (p p )} (1 + k 2 D 2 ) 

_ /2tt\ v l 

sin0> ( 


(t) • - = *-• 

\ A / u> uj 


(13) 


At Z = 300 Arm, the Debye length D e = V e /2u> 0 ~ 0.2 cm and because of the 
large wave length of the excited oscillations k 2 D 2 1 the second limitation 
of (13) is fulfilled when 


. Vi fum\ U! Vi 2 u 

sin 0 > — — = n ~ 10 2 — , 

il H \ c J ilfi c i Ih 


(14) 


and equation (12) becomes simpler and similar to the adequate equation in 
an isotropic plasma (0 ~ 0). Namely, 



■ - j%(pp )\ + 


(15) 


when 

Q2 

sin 2 0 < 1 , sin 0 <C 5.2 • 10 4 . (16) 

fi H 

In (14) and (15) the coefficient of refraction n = Vi = 1.5 • 10 5 cm/s and 
c = 3 ■ 10 10 cm/s. From (15) it follows that, at the altitude Z = 300 km, the 
value of < 5.2 • 10 -2 is very small. 
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Such ELF and VLF parametric oscillations were called by Aliev and Silin 
(see 1965 [7]) anisotropic sound waves. In our case in the ionosphere, they 
are ELF longitudinal waves propagating along the direction of the electric 
field Eq. Indeed, when a>2 lP 0// flo equation (15) becomes 


4 = \ {m) + = \ {0 < kE ”> 2 + u2v < = ( 17) 

where (see (4) and (5)) 

V. = 4-^tEo = ^4^ ■ Eo (18) 


x/2 


mu; 


is the velocity of this new parametric wave. 


The discovering of such waves in the ionosphere will be a very nice and 
an important contribution to the plasma physics. This can also be one ot 
the methods of diagnostic of the parameters of the plasma.The wave vector 
of these waves k = 27r/A may be estimated by these experiments. 


II. 2 Resonance oscillations in the frequency band 
swp = u>o, 0.25 < s < 2, s = 1 


Let us discuss here the behavior of the low frequency branch u ; 2 , p = (H p — > 
0), when the frequency of the electric field is not large and crosses the first 
resonance region iu p — su>o (s = 1). The formulas of the isotropic plasma 
are used for illustration of this dependence, particularly, because w 0 Wh- 
In this case, under the influence of the electric field the infinite spectra ol 
oscillation generated in the plasma are 


UJ = Wo 


1 5 = 00 

1 + 5 £ 

^ 5= — OO 


J 2 ,(Pr) 


ni 


(su> p + u> 0 ) : 


■) 


(19) 


and 



E J2 M 


(o> -F SUJpf ) 
(w + saj p ) 2 - oS 1 J ’ 


( 20 ) 


where in (19), (20) s = il, 2, — Taking into account that in (20) u> sui p , 
its denominator may be changed by |<su; p | 2 — an< ^ ^ ar enou gh from the 
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resonance regions where |su? p | 2 ~ u?q, the LF branch of waves (20) is estimated 
by the approximate equation 


7 rr 


UJ — Q,q [1 O r (/9p)] , $5 — , J r {pp)J -r{Pp) ? 

sin itr 


( 21 ) 


where r = ^ . In the resonance regions 


- 2 = ^A 2 i 

O 


1 ± 


1 


1 + 


(£)' 


( 22 ) 


where 


A s = 


, A p = 

32 


_V SU W 


W oJ 


1/3 


, A, < 1 (23) 


(Silin 1973 [9]). 


It is selfevident from equations (20 - 23) that by crossing the resonance 
regions ujq = u 2 p the frequency u> 2 becomes negative. I.e. that uj = 27 
and becomes the growth rate of the increasing disturbances of the plasma. 
Instabilities appear in the plasma when A s < 0, 0 < and > 1 (see (21)). 

W P 

In the resonance region A s < 1 and from (23) it follows that 



By using equations (21), the dependencies of and — on r = ^ sufficiently 
far from the resonance region, were calculated for Z = 300 km and p p — 1. 
The formulas (22) - (24) were used to calculate these dependencies around 
the resonance region u p = uo, 5 = 1. They are given on Fig. 2. Let us note 
here the following nice and important properties of this parametric branch 
of oscillations. 


When the frequency u p of the electric field Epo is comparable with the 
electron Lengmuir frequency, then in the resonance region, close to ^ ~ 1: 
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1 . The frequency fi p of the nonlinear parametric oscillations is no longer 
comparable with the ion Lengmuir frequency Q 0 (see Fig. 1). It reaches 
in our case the value fi p , ma;c = 9lO 0 by u 0 = 0.927u> p . Really it will be 
smaller in the collisional plasma (see II.3). 

2. A narrow region of an instability appears with a maximum of the 
growth rate 7 max = 44flo by u p = 0.972u;o- 7 = 0 when ^ = 1. 

3. A broad region of an instability appears with a narrow maximum 
Imax = 64flo by Up = 1.035u;o. Really, both the growth rates are 
smaller when the collision between the electrons and ions are taken 
into account (see next Section). 


To search and to find these wave branches and instabilities will be an 
important contribution to the TSS mission. 


II. 3 Excitation of oscillations close to the resonance region 


For excitation of the parametric oscillation, in particular, of the inter- 
esting for us here branches ~ — + o) , given on Figs. 1 and 2, the 

amplitude of the electric field is much smaller than when u p > u 0 ,u H close 
to resonance region. By the kinetic approximation it is determined by the 
condition 


0 < 


Au; = u p 


x/m+7®] < | vdi ■ vV + nj)- 


(25) 


Let us illustrate shortly this effect by the following. 

In the frequency band (25), the characteristic parameters of this process 
are described by the formulas (see [9]) 


r it = 


8^ e 

LJq 


■Dl 


k 2 = 


+ 2Au; 

TdTKT 


E l ~ 


kT 

D 2 VVeKie 


(26) 


where r^, cm is the amplitude of the oscillations of the electrons under the 
influence of electric field E 0 e‘ wt , v ei is the collision frequency between the elec- 
trons and ions and it is taken that N e nT e = ^/cT), D e — Di, and k 2 Dl <C 1. 
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When A u < 0, i.e. u p — yJujQ + fioj < 0 and the threshold values of Auj 

and Eq are determined by (26), the plasma becomes unstable and the process 
of exciting the parametric oscillations of the plasma begins. 


In the ionosphere, at the altitude Z ~ 300 km (see the Table in Section 

no 

i/ tj = 3 • 10 3 , NkT ~ 3.7 • 10“ 7 erg , D] ~ 0.2 cm 2 (27) 

and from (26) and (27) it follows that the threshold values of the character- 
istics of this process are 


1 « i _ £* = -15- 10 3 s~ l 

2(1- 3 k?Dl) ~ 2 

lO -5 cm 2 , E 0 = 3.83 • 1(T 5 CGSE = 1.1 Vfm 


and the maximum growth rate of these oscillation is estimated by 

[max - jg D 2 ~ 2 


Section III. Nonlinear heating of a magnetoplasma 

In the following parts of this section, we give the results of the theoretical 
study of the heating of a magnetoplasma under the action of the electric field 
Eoe*^ in the microscopic (hydrodynamic) approximation. They are used 
for some numerical calculations, mainly of the temperature dependencies of 
the electrons at different altitudes Z of the ionosphere and different angular 
frequencies a?, in particular, in connection with the orbit ( Z ~ 300 km) 
of the TSS-1 mission (Section IV.l). However, by the general statement of 
this problem, the microscopic theory has difficulties from the very beginning. 
The appropriate system of equations has not a stationary solution. This point 
becomes crucial at high altitudes of the ionosphere (Z > , 300 km) and, 

especially, in the magnetosphere and in a full ionized plasma, where the 
collision frequency between the ions and neutral particles v en ~ , = 0. In 
this case, the theory must be improved by introducing thermal losses in the 
adequate equations. Some of this effects are discussing shortly below (see 
III. 4). 
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III. 1 Statement of the problem. Microscopic theory 

Let us consider a magnetoplasma under the action of an alternative elec- 
tric field E = Eoe 1 ^, characterized at the starting point of this process t = 0, 
E 0 = 0 by the following. 

• it is an isothermal plasma, T e o = T,o = T n o - the temperatures of all 
the particles are equal, 

• it is a quasineutral plasma, N e0 = fV t0 - the densities of the electrons 
and ions are equal, and 

• it has one kind of neutral particles N n . The effective mass of these 

neutral particles and of the ions is equal to M n = Mi = , where 

N n = T,N S and N,, M s are the densities and masses of neutral particles 
of different kind. 

Then the temperatures and velocities of the electrons, ions, and neutral 
particles under the influence of E 0 e ,u;t can be calculated by the following 
systems of equations: 

_ _f_( Ve x B) - u ei {Y e - V.) - rUVe - V n ), 
m me 

% + -n x B ) + S-'-tv, - v.) - ^„(v f - v.), 

M Me M 

mN e u e n (V e - V n ) + - V n ), (30) 

and 

= _| c V e -E -6 ei v ei {T e -Ti)-6 en v en {T e -T no ) 
at 3 

^ = |eV,.E + 8 et u ei {T e -Ti)-6 in v in {Ti-T no ) 
at o 

^ n ~17 = N e 6 en u e n(T e -T no ) + NiuMT t -T no ) (31) 

at 

where 6 en , 6 ei , 8 in are, respectively, the energy lost by the electrons due to 
their chaotic collisions with the neutral particles and ions and energy lost 
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by the ions due to their collisions with the neutral particles. In the calcu- 
lation the values 6 en ~ £ ct ~ 2* 1CT 3 and Si n ~ 1 were used (see Gurevich [5]). 

In addition to equations (30) and (31), the following formulas of the 
collision frequencies must be used 



where i'en, o and are the initial values of the collision frequencies 

E 0 = 0, and the members ln[...] are the, so called, Coulomb logarithms. 

Thus, the systems of equations are interconnected and a selfconsistent, 
solution of (30 - 32) should be used to study in sufficient completeness the 
problem considered here. Namely this general statement is the crucial point, 
both for a full and correct understanding of the process of heating of a mag- 
netoplasma in the frames of the microscopic (hydrodynamic) approximation 
and also for a correct quantitative estimation of the temperatures. 


III. 2 Velocities V e j n * Approximation v — const 

By solution of the system of equations (34), the formulas of the velocities 
of the particles are visible when the dependence of the collision frequencies on 
the electric field is not taken into account - it is the approximation u = Const . 
The components of the electron velocities are the following: 


-eun 

v tx — Uy 

muj 


1 - i>f n - Cl 1 2 H + n ( u e { - j>e„) - i{nv ei + />,„) 

(1 - v e v in - uj 2 l ) + iv e -w|( 1 
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v ez 


— e 


E„ 


m u> 


i j(l + ii>in) 

(1 - VeVin) + Wt 

— q 2 h 

1 + (pen n ^*n u; +i„ t / ln 

} 


(1 - Ve^n -U>l) + 

2 

-ujf( 1 + IVin) 2 




rmo 


«(1 + ii>in) 

(1 - i'ei'in) + We 


(33) 


The components of the ion velocity are 


V = 

* IT 


Vi 


ly 


-ef l 


H 


mu) 


Ey i 




+ + P in ) - + ini', 


^cn, 


en u^-f tru/ t n 


2 + i u rn 




777.CJ 


j^(l - V e Vxn ~ &1) + tV e ~ Co 2 H { 1 + 0> ln ) 2 
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1 + i (p in - nu tn -^~) | 
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- U>//(1 + 2J> in ) 2 
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nv. 


eTl U/+tTU/| 


:)] 


(1 - VeVin) + iv e 


(34) 


and the velocity of the neutral particles V„, estimated by the velocities of 
the electrons and ions V e , V ; , is equal to 


V n = 


( v * + ;fc v .) 


(i + 


flUen n/lP e n ^ 


(35) 


The equations (33 - 35) are calculated in the coordinate system ( x,y,z ) 
where the magnetic field Ho = Ho(0,0,i/ 2 ) , i.e. it is parallel to the axis z , 
and the electric field E 0 lies in the plane ( yz ) i.e. E 0 = E o (0, E y , E z ), 

(0, E 0 sin 0, E q cos 0). Additionally, in (4) - (6) dimensionless collision fre- 
quencies 
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LJ 

U) 

u; 



(36) 


are used where fl//, lol — (Qh^h) 1 ^ 2 and ujh are, respectively, the angular 
ion, lower- hybrid and electron gyrofrequencies, /x = and n = 


III.3 Temperatures T e (E) and Tj(E) of the electrons and ions 


By calculating the temperatures, it is convenient and physically more 
acceptable to use the dimensionless values u/v instead of v ju> because u e = 
+ i> en is never equal to zero. Then, for example, the components of the 
velocities V etZ and Vi tZ (they are the same as in an isotropic plasma) are equal 
to 


K = V ez 
V, = v iz 


eEz 

mu e 

eE z 

~mV c 



^en 

V tn 



Uf 


2 


Ve »t n 




.(37) 


and their real parts are 

Re(V e ) 

Re(Vi) 


~e£ 0 !X e U; 2 - ~ Ve^n) 

m (u 2 - u e u , n ) 2 + v 2 lj 2 

— eEp VjJ_ - Ue n(u ; 2 - VeVjn) 

M (u 2 - u e u in ) 2 -\- i/]u> 2 


(38) 


Let us note here that by the approximation used in the earlier studies 
(see [5]), the following equations of the velocity V e of electrons 


^lU/t 


V, = 


m lJq + (jy € — iu) 2 


P / • \ . w //(EoH 0 ) w//(E 0 H 0 

E 0 (i / e - iu) + — H - 


Hl(v e - iu>) 


H o 


(39) 
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and of ,R e (V e ), when H 0 = 0, 


Re(V e ) = 


—eE o 

m 




UJ 2 + V 


2 ' 
e 


(40) 


were used. By comparison of (39) and (40) with (33) and (38), it follows 
that the equations of V e and certainly also of the temperatures, used in 
these studies, are effective only when 

UJ » v in , u 2 > v e u in , and u 2 > ui\ + v e v tn . (41) 


Visible, but rather complicated formulas of the temperatures of the elec- 
trons and ions T e (Eo,uj) and T{(Eo,uj) can be obtained only by the approxi- 
mation v = Const from the selfconsistent solution of the system of equations 
(30) and the shortened system of equation (31). Namely, when in the third 
equation of (35) ^ = 0. These formulas are as the following 


T. e 2 £ 0 2 

T n0 3m6 e is 2 T n0 


■*(-). 


(42) 


where 

$(...) = |i4i(j/) [cos 2 QR e (V ez ) + sin 2 0i2 e (K»)]} + 

4- {A 2 (v) [cos2u>tF eil (j/,u;) + sin 2u?t J F Ci2 (i/,u>)]} . (43) 


In formula (43), the members proportional to -jg • (R e { V ey , z}), and 
• Im{Vez] are omitted and the following notations are used: 


AM = 


A 2 {v) = 



(44) 


E e i(lS,Uj) 


F e2 {v,aj) 


— 2 ujB 2 -^ 3,2 B\ — 2uB2 

AI a + 4 uj 2 Al <2 + 4 to 2 

Az,\B 2 + 2uBi Az^B-j + 2ujBi 
Al tl + 4w 2 + Al 2 + iuj 2 ’ 


(45) 
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B\(u) = S e iy ei (l + -p^~] cos 2 QR e {V ez ] + sin 2 0/? e {K y }l + 

\ u e V / 

+ 2a; [cos 2 0/ m {V ez } -f sin 2 0/ TO { V ey } ] , 

B 2 {v) = s e u ei ^1 + -p^j cos 2 QI m {V ez } + sin 2 0/ m { V ey }] - 

- 2a; [cos 2 QR e {V ez } + sin 2 0/? e {V ey }] , (46) 

In addition, the values R e {V ey } , R e {V ez } , I m {V ey } and I m {V ez } are the real 
and imaginary parts of the terms in braces of the appropriate components of 
the velocity V e (see (33)). 

It is seen from equations (42, 43) that the temperature T e has two mem- 
bers. One of them estimates the stationary ( t — * oo) temperature, not de- 
pending on time. The establishment of the stationary temperature is illus- 
trated below by some results of calculations. The second member of (43) 
describes the harmonic periodical variation of the temperature; it is chang- 
ing in time by the double frequency of the electric field E ~ e' 2 *^ . In 
the following, we discuss only the stationary temperature. 


It is useful to introduce a characteristic field Eh (see [5]), which can be, in 
its own way, a measure of the intensity E l of the field needed for a remarkable 
heating of the magnetoplasma. When 0 = 0, i.e. in an isotropic plasma, 
from equations (42) and (43) it follows that 
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V e R e {V ez }] 1 • 


(50) 


where Ai(v) is given by (44) and 


A{v,u>) 


/ 1 \ 

2 7 

1 w 

m 

\ Win) 

(l ^ ' 

> + € 

\ ^e^in j 


Let us note here that by comparison of (42) and (49) we can see that the 
temperature Ti^Eq^u}) of the ions is smaller than the temperature T 6 (Eq,u>) 
of the electrons. However, really in the ionosphere this is remarkable only at 
the altitudes Z < (100 -> 200) km. At Z > 200 km, Ti(E 0 ,u) ~ T e (E 0 ,w). 


The time-dependent part of the ion temperature is determined by the 
same formulas as for the electrons (see (42, 43)), by replacing in the func- 
tions F(...) an £(...) (see (45) and (46)) the values V ey , V ez of the velocity 
components of electrons by the values Vi y and Vi z of the ions. The formulas 
(33) - (35) given above were used by some numerical calculations of the tem- 
peratures with the v = const approximation (see Section IV). Because of the 
diverse character of the problem studied here, let us briefly emphasize again 
some of its aspects. 


III.4 Some peculiarities of the problem. 

It is obvious that in the frame of the microscopic (hydrodynamic) approx- 
imation used in this study, the magnetoplasma should be a collisional media. 
However, if it would be, for example, a full ionized magnetoplasma, i.e. if it 
is characterized only by the constituent particles N e = N, (N n = 0) - by one 
kind of collisions v tx between the electrons and ions - then the used equations 
of the temperature would not have a stationary solution. The temperatures 
of the electrons would grow up continuously all the time of action of the 
electric field. It is because the loss of the energy of the electrons accelerated 
by the electric field is equal to the energy transmitted by them to the ions. 
Mathematically it means that the determinant A 2 of the first two equations 
of the system (2) becomes equal to zero. But this namely should happen, 
for example, when we consider these processes in the magnetoplasma at far 
distances from the Earth where u en — ► 0. 
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Within the limits of this study, the problem of the absence of a station- 
ary solution of the three equations of (31) exists from the very beginning if 
the equation of the derivative of the temperature of the neutral particles 
are taken into account. Mathematically it is because the determinant of the 
system (35) A 3 = 0. The temperature T n of the neutral particles is grow- 
ing up continuously under the action of the electric field, and can become 
even a source of heating of the electrons and ions. The crucial region of the 
ionosphere, where this process becomes highly noticeable is at the altitudes 
Z > 300 km (see below). The values T t and T t become larger in this region 
than they would be by the solution of the shortened system of equations (35). 

Thus it was important to learn in this study the applicability of the mi- 
croscopic approximation used in many studies and here by considering the 
behavior of all the temperatures T e , 71, and T n by the two approximations: 
v = const and i/(E). Especially it was important to search the role of the 
neutral particles. Let us note shortly here how the continuous growth of the 
temperatures of the plasma can be stopped. 

The heating of the magnetoplasma is accompanied and regulated, in ad- 
dition to the collisions n et1 u en and by many other processes. Their role 
should be of different degree. However, for a complete solution of this prob- 
lem, they should be taken into account. Then the general solution of this 
task will become stationary. One of these most important processes is the 
ionization of the neutral particles by the accelerated electrons under the ac- 
tion of the electric field E 0 e ,u,< . 

The velocity of the electrons becomes in the ionosphere, depending on 
the frequency u> and altitude Z, larger than the ionization potential E x of 
the neutral particles, namely, the thermal, chaotic velocity u e (E) plays the 
prime role by this process. It is considerably larger than the directed velocity 
V'e(E). For instance, in an isotropic plasma 

*(E = . *.WCE» . 2(WW ,E) (51) 

\ m ) vs 

(see (34) and (38)). The ionization potential of the atomic hydrogen H\ par- 
ticles (they are the main constituent at the high altitudes of the ionosphere) 
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is equal to £, = 13.54 eV. Therefore, the process of the ionization begins 
when Vje(E) = n e («) > 2.2 • 10 8 cm/s. It becomes sufficiently remarkable, i.e. 
that the additional density of the electrons N e ( E) can become even consid- 
erably larger than the regular density of the plasma N e when the amplitude 
of the electric field is rather small, even when E$ — (1 — 2) mV/m. This is 
seen from the following. 

The dependence of the cross-section cr e (Ei) of ionization of the atomic 
hydrogen on the energy of the electrons is given in Fig. 3 (Massy at al. 1969 
[12]). The value of cr e increases rapidly with E, up to a sharp maximum 
<7 e ,mar — (7 — 8)10 _17 cm 2 and then it is smoothly decreasing up to a eMat ~ 
2 • 10 -17 cm 2 . I.e. that at the altitude Z = 300 km, where the density of the 
atomic hydrogen N n ~ 3 • 10 9 , the production of the electrons 

/ e = <7 e • N n • V e (E) > (10 to 50) s” 1 (52) 

and the additional ionization is described by the recombination equation 

-j- = I e N — a e N 2 . (53) 

at 

It follows from (53), that growing up in time the electron concentration 
reaches a maximum by ^ = 0 equal to N = 10 9 , even if the coefficient of re- 
combination a e ~ 10~ 8 , cm 3 • $~V However, in this region of the ionosphere, 
a e < , < 10 -9 . Thus, under the action of the electric field, the surrounding 
plasma can become full ionized i.e. N e ~ 10 9 cm -3 . As a result, the collision 
frequency between the electrons and ions i/ ei (E) oc A r e (e)/T e (E) 3/2 becomes 
very large (see (3), and the T e (E) dependencies given in IV. 1). The growth of 
T e is stopped very quickly because the electron concentration is very quickly 
increasing. Indeed, in the point of inflection of the time dependence of the 
N(E, t), where = 0 , N = Thus, at that point ^ = Si 1S ver y 
large and the tangential of N(E,t) is almost vertical. Other losses of energy 
which can stop the growth of the temperature of the magnetoplasma are, for 
example, the following. 

The thermal emission of the constituent particles of the plasma i.e. the 
volume Stefan- Boltzman emission of the plasma can play a remarkable role. 
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The density of this emission is approximately equal to 


I p (ll>) = • n gr (u>)n 2 (uj), erg ■ s ■ cm 3 (54) 

In (54) n(tu) is a coefficient of refraction of the plasma and n gr = n -f u—J is 
the, so called, group velocity coefficient of refraction. 

The heating kinetic losses of the plasma should be considered and, in 
some cases, they can also be remarkable. I.e. the electron- electron and ion- 
ion collisions u ee and ua and possibly other kind of kinetic thermal losses 
which, in a sense, are similar to the Landau damping of e.m. waves in a 
plasma, should also be taken into account. 

To a certain extent, as it is seen below, we can omit in this work the 
discussion in detail of the noted above points and some other processes, ac- 
companying the heating of the magnetoplasma by an electric field. They are 
beyond the frame of this study. They are also not of the decisive importance 
at the altitudes of the ionosphere of the TSS Mission. 


Section IV. Numerical results. Conclusions 
Z = (10 2 to 10 3 ) km, F = (1 to 10 4 ) Hz 

In the following section, we give results of numerical calculations mainly 
of the temperature of the electrons and, for comparison, a few examples of the 
temperatures of the ions and neutral particles. The following model of the 
ionosphere (see Table) and the amplitude of the electric field Eq = 1 mV/m 
are used in this calculations. The altitude dependencies of the electron and 
ion densities and of the Lengmuir frequencies f 0 = u?o/27r and F 0 = Do/27r , 
given in the Table, correspond to the day-time conditions of the ionosphere. 
The altitude, frequency and angle 0 dependencies of the temperatures are 
presented in plots (Figs. 4 to 11). The frequencies F = uj/'lTt = 1 and 10 3 Hz 
and the angle between the magnetic and electric field 0 = 0 and the altitude 
Z = 300 km are selected by the calculations, because the orbit of the coming 
in the future TSS Satellite should pass at this altitude. 
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It is seen below, that the growth of the temperature with altitude in the 
frequency range F ~ (1 to 100) Hz is very large. It becomes about 10 2 to 10 3 
and more times larger than the initial temperature of the plasma already at 
Z ~ (200 - 300) km. Therefore, the frequency F = 1 Hz is used to illustrate 
more dramatically the region where it is necessary to stop the temperature 
growth by including in the theory additional losses of energy (see above III. 4). 
The zone of ionization of the plasma by the accelerated electrons is therefore 
marked on the plots by a thin line and vertical arrows. The frequency band 
F = (1 to 10 4 ) Hz, that we have used, is much smaller, but also becomes 
somewhat larger than the gyrofrequency Fh — ^ of the ions in the altitude 
range discussed. It is also smaller and much smaller than the low hybrid fre- 
quency Fi = In addition, the value = (v e i + Vi n )vi n is much larger 

than this frequencies at low altitudes and becomes comparable and also much 
larger than at high altitudes. Therefore, the formulas used in the literature 
are ineligible for these calculations (see page 18 and (41)). The used angle 
0 = 0 corresponds to the isotropic plasma and characterize the behavior ol 
the maximum values of the temperature. However, the evaluation of these 
dependencies on different values of 0 are also given and the influence of the 
magnetic field H 0 is shown by these data. 

The results of the numerical calculations mainly correspond to the ap- 
proximation v = const , i.e. that the shorted system of (31) is used and it 
is assumed that = 0, r n = r„„. However, results of the selfconsistent 
solution of the full systems of (30) to (32) are also given by some examples 
and the role of the neutral particles is discussed. 

Let us now describe in more detail the numerical results by the plots, 
presented here . 
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IV. 1 Altitude dependencies of temperatures. 


a. 


dTi 

dt 


= 0 


b dla 
D ' dt 


# 0 . 


a.^“- = 0. First of all let us consider the altitude dependencies of the 
ratio of the temperatures T e (Eo,u) and Ti(E 0 ,u>) to the initial temperature 
of the plasma T n0 , shown on Fig. 4 for F = (1 and 10 3 Hz , calculated by the 
shorted system of equations (31) and by the two approximations: u = const 
and v — v(T) = v(E 0 ,u>). How the collision frequencies i/(E,<) reaches their 
stationary values at the altitude Z = 10 3 km and F — 1 Hz , is also shown 
for sake of illustration, on Fig. 5. 


The following is seen by examination of this figure. 

1. The temperatures of the electrons and ions are growing up quickly 
on both frequencies. They become equal in the altitude region Z ~ 
(200 - 250) km. 

2. The temperatures, when F = 10 3 Hz, have a maximum at Z ~ 200 km 
by the approximation u — v(E 0 ,uj). The same happens by approxima- 
tion v = const when F = 1 Hz. 

3. At Z > 300 km the temperatures are smaller by the approximation 
v = n(E 0 ,u>) than by the approximation n — const. 

4. The process of ionization by the accelerated electrons is acting almost 
in all the altitude region when F = 1 Hz and does not play a role , 
when F = 10 3 Hz, especially at Z > 300 km. Thus, the process of the 
ionization should be stopped when F = 1 Hz already at the altitude 
Z = (120 - 150) km. 

The properties noted above of the nonlinear heating process of the iono- 
sphere are typical. Namely, the behavior of T e and 7i for F = 1 Hz is 
typical for the frequency band F = (1 to 10 2 ) Hz and at frequencies a lit- 
tle higher. The behavior for F = 10 3 Hz is typical for the frequency band 
F > (10 2 — 10 3 ) Hz up to the value F = 10 4 Hz used here. This is seen from 
the following dependencies of T e (Eo,u>)/T n o on Z: 

• in Figure 6, by the curves calculated by the approximation v = const 
for F = (1, 10 2 , 10 3 andlO 4 ) Hz, 
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• in Figure 7, by the curves for the angles 0 = (0, 30, 60, 75, 85, and 90) 
degrees also calculated by the approximation v — const , 

• in Figure 8, by the curves, calculated for F = 10 3 Hz and the angles 
0 = (0 and 75) degrees by the v = u(E 0 ,uj) approximation. 

The characteristics of these altitude dependencies of the temperature arc 
in agreement with the noted above points 1 to 4. 


^ 0. Some results of the selfconsistent solution of the full systems 
of equations (30)- (32) are presented in Figure 9. Not to overload this figure, 
only the T e and T n dependencies are given in it. The main new important 
characteristics of these dependencies in addition to the given above in a. of 
this section are the following 

1. The temperature T n of the neutral particles is growing up very quickly 
at altitudes Z > 200 km. It approaches the temperature T e of the 
electrons at Z > (400 — 500) km and becomes close to T e similarly to 
the temperature of ions, T t . 

2. At the altitudes Z > 400 km, the temperatures T e are larger on both 
the frequencies by the approximation — ^ ^ 0 than by the approxima- 
tion a* = 0. 

3. At least at Z > 300 km, the heating of the neutral particles becomes 
a source of the heating of the electrons and ions. The growth of T n 
should be stopped, similarly to the case = 0 at these altitudes (see 
III. 4). 

4. In general, the heating of the neutral particles is an important part in 
many cases of the process of heating of the plasma. 

How the temperature T c (Eq, t) reaches the stationary temperature at the 
altitude Z = 200 km is shown on Fig. 10 for F = 1 Hz and 0 = 0. It is 
seen how slow this process is going on when ^ / 0, in comparing with the 
case to = 0. 
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IV. 2 Frequency dependencies of the temperature. Z = 300 km. 

For the interested to us altitude Z ~ 300 km, the frequency dependencies 
of T e are given on Fig. 11. They were calculated for the both approximations 
and E 0 = 1 mV/m. In a sense, this is a a transition region. The influence 
of the heating of the neutral particles is still not too large in this region (see 
Fig. 9). 

The behavior of the frequency dependencies of T e is presented on Fig. 11, 
for different angles 0, are in agreement with the given above conclusions in 
the sections IV. 1. a and IV.l.b. The heating of all kind of particles is very 
high in this region. Even by amplitudes of the field E 0 ~ (10~ 2 — 10' - l ) V/m, 
the temperatures of the electrons and ions can reach rapidly, at least, values 
T e (E 0 ,u) > (10 - 10 2 )T nO in a broad range of angles 0 between the electric 
field E 0 and the Earth’s magnetic field H 0 . Temperatures T e {E 0 ,u) > 10 2 r n0 
should, at least, be stopped by the ionization and recombination and other 
processes, accompanying the heating of the plasma by the electric field E = 
E 0 e iwt . This does not mean that the temperatures of the electrons and ions 
can not become larger than the critical temperature T e< (,) of ionization, let us 
say, of the atomic hydrogen#!, which is equal to 

T e>(l) = 13.5 eV = 1.16 • 10 4 • 13.5 = 1.57 • 10 5 , I<°. (55) 

However, to answer this question, theoretical calculation must be done taking 
into account losses of energy, in particular, noted in Section III. 4. 


Summary 

Some nonlinear parametric and heating effects in a magnetoplasma were 
considered in this work, particularly, in connection with the forthcoming 
Tethered Satellite Mission (TSM). For this aim, the theory of nonlinear heat- 
ing of a collisional magnetoplasma, was extended, taking into account the 
heating of all the kinds of particles and all the kind of collisions between 
them and also the interconnection of these processes. It seems to the author, 
that at first the generating of parametric branches of oscillations due to the 
decay instabilities are discussed first in this work in connection with iono- 
sphere. This was done on the basis of the known theories. 
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The numerical results describe in detail the altitude, frequency etc. tem- 
perature dependencies. Conclusions, as they follow by describing of different 
aspects of the discussed problems, are given. 

These phenomena can play a great role in the plasma surrounding the 
TSS System. It will be an important contribution both to plasma physics 
field and for diagnostic of many properties of the plasma to search and to 
investigate them in this unique experiment. By the theoretical treatment of 
the related experimental data, the wave numbers and other local parame- 
ters of different branches of waves can be estimated. The processes of the 
growth of generation of different branches of e.m. oscillations and waves and 
of the artificial ionization of the plasma, and many other phenomena can be 
learned as they are developing in the course of this experiment. Despite of 
the fact, that the TS System will not have special artificial sources of alter- 
nating electric fields, the electromagnetic oscillations, generated in the TM 
Cloud by other kinds of instabilities can become such sources. Possibly, they 
will produce heating and parametric effects. 

The great importance, for many reasons, to search for these effects in 
the ionosphere and magnetosphere, makes it advisable, that generators of 
alternating electric fields in the ELF, VLF and HF bands should be added 
to the TSS System in the future TSS Missions. 
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Figure 1. 

Resonance branches of a collisionless magnetoplasma 
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Figure 3. Experimental and theoretical dependencies of the cross-section 
a e on the potential E ly eV of ionization of the atomic hydrogen H j by 
accelareted electrons (Massey et al. 1969 [12]). 






Figure 6. The same as on Fig. 4. Approximation V = const ,0 = 0. 
1. F = 1 Hz , 2. F = 10 2 Hz, 3. F = 10 3 Hz 
4. -* F = 10 4 Hz 
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Figure 7. The same as on Fig. 4. Approximation V = const . F = 1 //;. 
1 . , 2 . , 3 . , 4 . , 5 . , and 6. -» 6 = 0 , 30 , 60 . To 
S5 . 90) degrees between the vectors E,, and H 0 . 




Figure 8. The same as on Fig. 4 for F = 10 3 Hz. 

Approximation V = const : lo— +0 = O:2o— ^0 — ” ;) 
Approximation V = y(£o^) : 0 — 0 » 0 = i ’ ) 



Figure 9. Altitude dependencies of the ratios of the temperatures^ 
Approximation V = V(Eq, uj) ,0 = 0. 

Zff 1 : F = 1 Hz, 1. -» Zfr = 0 , l T - ^ * o 

F = 10 3 Hz , 2. - = 0 , 2 T - 1 0 

: F = 1 ffz, 3. - ^ = 0 , 3 T - ^ # 0 


t, sec : a x 10 2 , b x 10 3 


Figure 10. The establishment of the temperature T ( {E 0 .~) in time. 

F = 1 Hz , 0 = 0, Z = 200 km. Approximation V = U{ E 

a. 
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Figure 11. Frequency dependencies of the ratios Te ^~ ^ ■ Z - 300 m ■ 

Approximation V = const : 1. -* 0 = 0 ; 2. — >0 = 75° . 3. — <• 0 = S)° 
Approximation V = V(Eq,^) ■ 1e — * © = 0 ; 2g — ► 0 = 75° ; 3^ — » 0 = 8 > 





